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Abstract 



The usefulness in control theory of the geometric theory of motion 
on Lie groups and homogeneous spaces will be shown. We quickly 
review some recent results concerning two methods to deal with these 
systems, namely, a generalization of the method proposed by Wei 
and Norman for linear systems, and a reduction procedure. This last 
method allows us to reduce the equation on a Lie group G to that on 
a subgroup H, provided a particular solution of an associated problem 
in G/H is known. These methods are shown to be very appropriate 
to deal with control systems on Lie groups and homogeneous spaces, 
through the specific examples of the planar rigid body with two oscil- 
lators and the front-wheel driven kinematic car. 



Key words: Drift-free control systems, Wei-Norman method, motion in Lie 
groups and homogeneous spaces, reduction. 
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1 Introduction 



Mechanical systems whose configuration space is a Lie group are known to 
be of great relevance both in mechanics and control theory. The simplest 
example is that of a point particle moving freely on M 3 , for which the con- 
figuration space can be identified with the Abelian group of translations in 
a three-dimensional space. The motions of the rigid body are described by 
the Euclidean group in three dimensions, E 3 = T 3 £70(3, R), and for the 
rigid body with a fixed point the configuration space is 5*0(3, M). However, 
there are also other systems which appear often in control theory which can 
be formulated in homogeneous spaces and reduced to problems of motion on 
Lie groups. The techniques used here have been developed in (21 El HI IE] , and 
firstly applied in control theory in GH| • The paper is organized as follows. 
In the next section we review briefly the geometric theory of motion on Lie 
groups and homogeneous spaces. In Section 3 we recall two basic methods for 
dealing with such systems, and in Section 4 we illustrate the theory through 
its application on two control systems, namely, the planar rigid body with 
two oscillators and the front-wheel driven kinematic car. 



2 Motion on Lie groups and homogeneous 
spaces 

Given an r-dimensional connected Lie group G, the set of curves 7 : K — > G, 
t 1 — > g{t), can be endowed with a group structure by means of the point- 
wise composition law, 72 * 71 : t ^ 92{t)gi{t), for all { 6 1, Moreover, 
G acts by right and left translations on itself, and if g(t) defines a curve 7 
starting from g(0) = e, then, for each g G G, the maps 7 * 70 : t 1— > g(t) g 
and 70 * 7 : t 1— > gog(t) are two curves starting from g , which are called, 
respectively, the right and left translation of 7 by g . 

In addition, the map 7 : K — > G given by t 1 — ^ g(t) defines a vector 
field along the curve g(t), and then, using the right translation by we 
obtain a curve in T e G which can be written as 

r 

g(t)g-\t) = -J2b a (t)a a , (1) 
0=1 

where {a 1; . . . , a r } is a basis of the tangent space T e G at the neutral element 
e G G. The minus sign on the right hand side is a matter of convention. 
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The left hand side should be understood as R g -iM* g M(g(t)), although for 
the particular case in which G C GL(n,M), for certain n, it reduces to the 
above expression. 

A first important point is that the curve g(t) can be recovered as the only 
solution of equation starting from e G G. Any other solution is obtained 
from the previous g(t) by right translation, and, in this sense, (P) is right- 
invariant. In other words, g(t) is the integral curve starting from e G G of 
the t-dependent vector field in G 



where is the right-invariant vector field in G whose value in e G G is a a , 
Xa (g) = R g *e(o, a ). Similarly, X^(g) = L gJre (a a ) denotes the left-invariant 
vector field on G determined by a a . The left-invariant vector fields in G 
close on a finite-dimensional Lie subalgebra q of %(G), which is called the Lie 
algebra of G. The right-invariant vector fields close on a finite-dimensional 
Lie algebra opposite to g. 

Let us consider now a transitive action $ : G x M — > M of G on a 
manifold M (which is then called a homogeneous space of G). Each choice 
of a point xo G M allows us to identify M with a space of left cosets, i.e., 
M = G/H, where H is the stability subgroup of the point x with respect 
to Different choices for x lead to conjugate subgroups. We recall that G 
can be regarded as the total space of the principal bundle (G,r,G/H) over 
G/H, where t:Gh G/H is the canonical projection. 

The important point now is [3] that the right-invariant vector fields 
are r-projectable, the projections being the fundamental vector fields in M 
associated to the natural left action of G on M = G/H, identified with $. 
More explicitly, we have X a (x) = X aa (x) = & x * e (— a a ), where x = gH, and 
T* g Xa(g) = —X a (gH). Consequently, the vector field (0) projects onto the 
t-dependent vector field in M 



r 




(2) 



r 




(3) 



giving rise in this way to a system of differential equations 



r 




I = 



1, . . . , n = dimM . 



(4) 



a=l 
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This relation between both systems tell us that the solution of (j3J) starting 
from the point x G M, arbitrary but fixed, is given by x(t) = $(g(t),x ), 
where g(t) is the solution of (0) starting from the identity. Note that the 
vector fields arising in (j3J) close on the Lie algebra g. 

Conversely, if a system of type (jU) is defined by complete vector fields 
closing on a finite dimensional Lie algebra g, then it can be regarded as a 
system of the described type, taking a Lie group G such that the vector fields 
are the fundamental vector fields of its action on the manifold M given by 
the integration of their flows. 

In this sense equation (JTJ has a universal character, since the knowledge of 
its solution starting from the identity is enough to write the general solution 
of each associated system of type (J3J) on each homogeneous space for G. 
Unfortunately, finding the desired particular solution of (0) may be a difficult 
problem. We briefly recall here two methods for dealing with such equations. 
The first one is based on a generalization [JJ |S] of the method proposed by 
Wei and Norman for linear systems [TTl I18j . The second one is a reduction 
procedure allowing us to reduce the problem to a similar one in a subgroup, 
when a solution of the associated problem in the corresponding homogeneous 
space is known E3 HI IE] • 

3 The generalized Wei— Norman method and 
the reduction procedure 

The main idea of the generalization of the Wei-Norman method consists on 
writing the curve g(t) solution of (JTJ starting from the identity, in terms 
of the second kind canonical coordinates with respect to a basis of the Lie 
algebra g, {ai, . . . , a r }, for all t, i.e., we write g{t) = Yl a =i ex P( — v a {t)a a ) . 
Then, the differential equation (0) transforms into a system of differential 
equations for the v a (t), and we have to find the solution determined by the 
initial conditions v a (0) = for all a = 1, . . . , r. A simple calculation shows 
|H1 that 



R g (t)-i * g (t)(g(t)) 



a=l \/3<a 
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and then, substituting into equation (JTJ) we obtain the fundamental expres- 
sion of the Wei-Norman method 

r / \ r 

y^t>q I TT exp(-f /3 (t) ad(a /3 )) J a a = ^ b a (t)a a , (5) 

a=l \/3<a / a=l 

with f a (0) = 0, a = 1, . . . , r. The resulting system of differential equa- 
tions for the functions v a (t) is integrable by quadratures if the Lie algebra is 
solvable [HI EE], and in particular, for nilpotent Lie algebras. 

On the other hand, given an equation like (Q) on a Lie group G, it may 
happen that the only non- vanishing coefficients b a (t) are those corresponding 
to a subalgebra f) of g. In that case, the equation reduces to a simpler equation 
on a subgroup, involving less coordinates. 

Now, the important result is that the general situation can be reduced 
to this simpler one under certain conditions. In fact, we can show that the 
problem of finding the solution of starting at e G G can be reduced to 
the one of solving a similar equation in a subgroup H, provided that one 
particular solution gi(t) for the system defined by the left action $ of G on 
the homogeneous space M = G/H is given. The result is as follows j3]: 

Theorem 1 Every integral curve of the time- dependent vector field HJ) in 
the group G can be written in the form g(t) = gi(t)h(t), where gi(t) is a 
curve projecting onto a solution gi{t) of the system of type associated to 
the left action $ on the homogeneous space G/H, and h(t) is a solution of 
an equation of type (QJj but for the subgroup H , given explicitly by 

(hh~ l ){t) = -Ad(g?(t)) (^Ut^ + ^g^ 1 )^ e T e H . (6) 

This result allows us to understand from a group theoretical point of view, 
as an example, the classical integration theorems for the Riccati equation 
involving one, two or three of its particular solutions [7j, see also jUE]- 

4 Some illustrative examples in control the- 
ory 

We will apply the previous techniques to problems appearing in the context 
of control theory, after a very brief description of such kind of problems. 
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Roughly speaking, a control system is a dynamical system which depends 
on a set of control functions, u a , which model magnitudes to be varied some- 
how externally to the system, in order to modify its behaviour according to 
a specific purpose. For example, it could be desired to reach one point from 
another one of the configuration space, or to minimize some cost functional 
along the evolution of the system with prescribed initial and final conditions. 

From a geometric point of view, we can describe the framework of a 
control system as a bundle B (usually a trivial vector bundle) on the state 
space manifold M, with projection ttb '■ B — > M. The control dynamical 
system corresponds then to the integral curves of a vector field along the 
projection tt b , which in local coordinates, (x\u a ), reads 

x i = X\x, u), x E M , u E B x = TTg^x) . (7) 

A special class is those of the control systems affine in the controls: 

r 

x i = X l {x) + J2 U *K(%) , 2 = l,...,n = dimM. (8) 

a=l 

The vector field X l Q d/dx l G X(M) is called the drift of the system. In 
principle, the control functions u a are supposed to depend on time (in the 
control theory terminology the system is then operating in open loop). 

We are now interested in those which are drift-less, i.e., X = 0, since they 
become a system of type as soon as the set of vector fields {X l a d / dx 1 } 
closes on a finite-dimensional Lie algebra. 

It may also happen that starting with a control system that is of the form 
(J3J), but where the vector fields {X l a d / dx 1 } do not close on a Lie algebra, an 
appropriate feedback transformation u a (x,t) = Y^p=i fa0{ x ) v fi(t) could lead 
to a new system written as 

r 

= ^2v a (t)Y*(x) , i = 1,..., n = dimM, (9) 

a=l 

where Y a (x) = Ym=i fpot{.x)Xp(x) do close a finite-dimensional Lie algebra. 

We will illustrate by means of simple examples how systems of type 
both directly or as a consequence of a feedback transformation, appear in 
control theory. In addition, with the general theory developed in previous 
sections, we will show how it is possible to solve, reduce, or relate systems of 
this kind formulated on different homogeneous spaces. 
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4.1 Planar rigid body with two oscillators 

This example comes from the consideration of the optimal control problem of 
a planar rigid body with two oscillators JH]- The control system of interest 
has the state space R 2 x S 1 , with coordinates (x\, x 2 , 6) 

xx = b^t), x 2 = b 2 (t), 9 = x\b 2 {t) - x 2 2 b x {t) , (10) 

where bi(t) and b 2 (t) are the control functions. This system is similar to the 
celebrated Brockett nonholonomic integrator system P, but where the third 
equation is quadratic in the coordinates instead of linear, and the meaning 
of the third coordinate is now an angle. 

The solutions of the system ()1U|) are the integral curves of the time- 
dependent vector field b\{t) X x + b 2 (t) X 2 , with 

The Lie brackets 

X 3 = [X,, X 2 ] = 2( Xl + x 2 )—, X 4 = [X u X 3 ] = 2— , 

jointly with X\, X 2 , make up a linearly independent set in points with X\ ^ 
—x 2 , and the set {X l; X 2 , X 4 } spans the tangent space at every point of 
M 2 x S 1 . According to Chow's theorem every two such points can be 
joined by appropriate piecewise constant controls &i(t) and b 2 (t), therefore 
the system is controllable. In addition, the set {X\, X 2 , X 3 , X{\ closes on 
the nilpotent Lie algebra defined by the non-vanishing Lie brackets 

[X 1: X 2 \ = X 3 , [X h X 3 ] = X 4 , [X 2 , X 3 ] = X 4 , (12) 

isomorphic to a nilpotent Lie algebra, denoted as g 4 , which can be regarded 
as a central extension of the Heisenberg Lie algebra f)(3) by R. In fact, taking 
the basis {ai, a 2 , a 3 , a 4 } of g 4 for which the non-vanishing Lie products are 

[ax, a 2 ] = a 3 , [a l5 a 3 ] = a 4 , [a 2 , a 3 ] = a 4 , (13) 

then the center 3 of the algebra is generated by {a 4 }, and the factor Lie 
algebra g 4 /$ is isomorphic to the Heisenberg Lie algebra f)(3). 
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Let G 4 be the connected and simply connected nilpotent Lie group such 
that its Lie algebra is the previous g 4 . The right-invariant system of type (0) 
on G4 corresponding to (fTUjl is 

Rg(t)-i* g (t)(g(t)) = -b\{t)ai - b 2 (t)a 2 . (14) 

Let us solve it by the Wei-Norman method. We write the solution of (|14j) 
starting from the identity as the product of exponentials 

g(t) = exp(-t>i(t)ai) exp(-t> 2 (t)a 2 ) exp(-w 3 (t)a 3 ) exp(-v 4 (t)a 4 ) , (15) 

and using the expression of the adjoint representation of g 4 , and its exponen- 
tiation, and applying then ©, we find the system 

u 1 = 6i(t), v 2 = b 2 (t), v 3 = b 2 (t)v 1 , v 4 = b 2 {t)v 1 (v 1 /2 + v 2 ) , (16) 

with initial conditions v\(0) = v 2 (0) = 1*3(0) = 1*4(0) = 0. The solution is 
easily found by quadratures: if we denote Bi{t) = J^b^s) ds, i — 1, 2, then, 



Vl (t) = B 1 (t) , v 2 (t) = B 2 (t) , v 3 (t) = [ b 2 { S )B 1 ( 8 ) ds . 

Jo 




Other orderings in the factorization (fT5j) are possible, with similar results. 

Now, we can find the expressions of the action $ of G4 on the con- 
figuration manifold R 2 x S 1 such that Xi be the infinitesimal generator 
associated to a,i for each i 6 {1, . . . , 4}, and of the composition law of 
G4. We will use canonical coordinates of the second kind in G 4 defined 
by g = exp(aai) exp(6a 2 ) exp(ca 3 ) exp(da^), and then, integrating and com- 
posing accordingly the flows of the vector fields X iy the action reads $ : 
G 4 x (R 2 x S 1 ) —>R 2 x S\ 

$((a, b, c, d), (xi, x 2 , 0)) = (xi - a, x 2 - b, 

9 + axj- bx\ -2(ab + c)x 2 - 2cx x + ab 2 - 2d) , (18) 

and the composition law of G4 reads 

(a, b, c, d)(a', b', c', d') = (a + a', b + b', c + c' — bo! , 

d + d' - c(a' + b') + ba\b + 2b' + a!)/2) . 
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The neutral element is represented by (0, 0, 0, 0) in these coordinates. 

The general solution of (fTU|) can be calculated by means of the solution 
of the Wei-Norman system (fTBj) as 

$((-t>i, -v 2 , ~V 3 , -V4), (X10, X 20 , 9q)) = (X10 + Vx, x 20 + v 2 , 

9 - v x x\ Q + v 2 x 2 w - 2{v x v 2 - vz)x 20 + 2v 3 x w - v x v\ + 2v 4 ) , 

where V\ = t>i(t), v 2 = v 2 (t), v 3 = v 3 (t) and f 4 = v^it) are given by (|T7|) . 
(xio, x 20 , 9 ) e I 2 x S 1 are the initial conditions and $ is given by (fTS|) . 

We will show the way the previous reduction theorem applies to the con- 
trol system (fTUj) . Several possibilities of reduction exist, as many as (non- 
equivalent) subgroups of G 4 . Specifically, we will show how the system (fPIj) 
can be reduced to a control system of Brockett type plus a system on the 
real line, performing the reduction with respect to the center of G4. 

For a better illustration we parametrize now G4 by canonical coordinates 
of first kind defined by g — exp(aai + ba 2 + ca 3 + da^), being then the 
composition law 

(a, b, c, d)(a', b' , c', d') = (a + a, b + b' , c + c + (ab' — ba')/2, 
d + d'+ {ad - ca')/2 + (be - cb')/2 + (a&' - ba')(a -a' + b- b')/l2) . 



Thus, the adjoint representation of the group reads 



/ 



Ad(a, b, c, d) 



1 





1 

a 





1 



\ 






(19) 



\-|(a + 6)-c |(a + 6)-c a + 6 1 / 



If g{t) = (a(t), b(t), c(t), d(t)) is a curve in G 4 expressed in the previous 
coordinates, we obtain 



/ 



a 
b 



\ 



c — \ (bd — ab) 



\d - Uab + b 2 + 3c)d + Ua 2 + ab - 3c)b + \(a + b)c J 



To perform the reduction we select the subgroup H to be the center of 
G 4 , which is generated by a 4 . The relevant factorization is (a, b, c, d) = 
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(a, b, c, 0)(0, 0, 0, d). Therefore, we can describe the homogeneous space 
M = G A /H by means of the projection 

t : G 4 — ► GJH 
(a, b, c, d) i — >■ (a, 6, c) , 

associated to the previous factorization. We take coordinates (yi, y 2 , y 3 ) in 
M. Thus, the left action of G A on such a homogeneous space reads 

$ : G A x M — ► M 
((a, 6, c, d), (y x , y 2 , y 3 )) 1 — ► (3/1 + a, y 2 + b, y 3 + c + {ay 2 - by x )/2) . 

The corresponding fundamental vector fields are 

yH - _ 9 _ 2/ 2 9 yH _ _ 9 ,2/1 9 vH _ _ 9 Y H _ n 
1_ flfc 2^3' * 2 ~ dy 2 + 2dy 3 > As " ^ ~°' 

which span the tangent space at each point of M, and, in addition, satisfy 
the Lie brackets [Xf , Xf ] = X* , [Xf , Xf ] = Xf = and [Xf , Xf ] = 
Xf = 0. 

Now, let the desired solution of be factorized as the product 

gi(t)h(t) = ( Vl (t), y 2 (t), y 3 (t), 0)(0, 0, 0, d(t)) , 

where gi(t) projects onto the solution r(gi(t)) = (yi(t), y 2 (t), y 3 (t)), with 
initial conditions (0, 0, 0), of the Lie system on the homogeneous space M 
associated to ([14)1. 

£i = -6i(t), y 2 = -b 2 (t), y 3 = -{b 2 {t) yi -b x {t)y 2 ). (20) 

Then, we reduce the problem to a Lie system in the subgroup H for h(t) = 
(0, 0, 0, d(t)), with /i(0) = e, i.e., c?(0) = 0. The expression of this last 
system is given by Theorem ^ i.e., 

d = —((2/1 + 2/2) - &22/1) - Qysih + b 2 )) , 

which is a system on H = R, solvable by one quadrature. The system 
(j20|) is of Brockett type [l] (indeed they are related by the simple change of 
coordinates x = —y±, y = —y 2 and z = —2y 3 ), and therefore we obtain two 
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interesting results. Firstly, that solving a system of type (jTUjl can be reduced 
to solving first a system of Brockett type and then to solving a Lie system in 
M, which is immediate. Secondly, that the Brockett system can be regarded 
as a Lie system on H(3) written, moreover, in terms of canonical coordinates 
of first kind of such a group. 

As an interesting open problem, it remains to investigate the interrelations 
the corresponding optimal control problems might have with respect to this 
reduction. 



4.2 Front- wheel driven kinematic car 

The front-wheel driven kinematic car has been considered by a number of 
authors, mainly with regard to the nonholonomic motion planning problem, 
and as such is made nilpotent by a state space feedback transformation 12, 
El El El E] • It is a simple model of a car with front and rear wheels. The 
distance between the rear and front axles is I, which we will take as 1 for 
simplicity. 

The configuration of the car is determined by the Cartesian coordinates 
(x, y) of the rear wheels, the angle of the car body 9 with respect to the 
horizontal coordinate axis, and the steering front wheel angle (ft £ I — 
(— 7r/2, 7r/2) relative to the car body. The configuration space is therefore 
R 2 x S 1 x 7, with coordinates (x, y, 9, <fi). The controls of the system are 
the velocity of the rear (or sometimes front) wheels and the turning speed of 
the front wheels. For a schematic picture of the system, see, e.g., [T3} ITE]. 

The control system for this model can be written as El (compare 
with [HI Eq. (13.7)]) 

x = ci(t) , y = c\{t) taxi 6 , <fi = c 2 (t) , 6 = c\(t) tan0sec# . (21) 

Note that this system is defined for angles 9 with cos# 7^ 0. We therefore 
restrict 9 G I as well. The solutions of (|21j) are the integral curves of the 
time-dependent vector field C\{t)Yi + c 2 (t)Y" 2 , where 

Y x = ^- + tan^ + tan0sec^, Y 2 = ^-. (22) 
ox oy o9 o<p 

Taking the Lie brackets 

Y 3 = [Yi, Y 2 \ = -sec9 sec 2 <f>— , Y A = [Y 1: Y 3 ] = sec 2 9 sec 2 0— , 
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we see that {Yi, Y 2 , Y 3 , Y4} generate the full tangent space at points of the 
(restricted) configuration space M? x I x I, thus the system is controllable 
there. However, (j2"Tj) is not a system of type (j3J), since the iterated Lie 
brackets 

[Y 2 , [Y 2 , . . . [Y 2 , Y x ] ■ ■ ■ ]] or [Y u [Y u . . . [Y u Y 2 }--- ]] 

generate at each step vector fields linearly independent from those obtained 
at the previous stage, therefore they do not close a finite-dimensional Lie 
algebra. 

Notwithstanding, it can be transformed into a nilpotent system of type 
(jlj). In [121 EH HE] it is proposed the following state space feedback trans- 
formation (it seems that in [TH ITH] there are some minor misprints, for their 
expressions do not do the work) 

ci(t) = &i(t) c 2 (t) = -3 sin 2 sec 2 ^ sin ^ 61 (t) + cos 3 # cos 2 </> & 2 (t) , (23) 

and then the change of coordinates 

xi = x , x 2 = sec 3 9 tan <fi , x 3 = tan 9 , X4 = y , (24) 

with inverse 

x 2 



x = X\ , y = £4 , 9 = arctan x 3 , = — arctan 



(1 + xlf' 2 



which transforms (J21|) into the control system in R 4 with coordinates (x 1; x 2 , £3, £4) 
given by 

x 1 = b 1 (t), x 2 = b 2 (t), x 3 = bi(t)x 2 , x 4 = h(t)x 3 , (25) 

where the control functions are now bi(t) and b 2 (t). In other words, we take 
the new input vector fields 

X x = Yi - 3 sin 2 <p sec 2 9 sin 9 Y 2 , X 2 = cos 3 9 cos 2 <p Y 2 . 

The system is usually said to be in chained form, see [HJ [TH] . Let 
us show that it is a system of type (@J). The new input vector fields read, in 
the coordinates (xi, x 2 , x 3 , £4), as 

d d d ^ d , . 

X 1 = — + x 2 —+x 3 —, X 2 = —. 26 
OXi OX 3 OX4 ox 2 

13 



The Lie brackets 

d d 

X 3 = [Xi, X 2 ] = 5 — , X A = Xl, X 3 ] = — — , 
0x3 0x4 

are linearly independent from X\ and X2, and {X\, X 2 , X 3 , X 4 } generate the 
full tangent space at every point of the configuration space M 4 , so the system 
is controllable. On the other hand, the same set closes on the nilpotent Lie 
algebra defined by the non-vanishing Lie brackets 

[X 1 ,X 2 ]=X 3 , [X 1 ,X 3 ]=X A . (27) 

This Lie algebra is isomorphic to a four dimensional nilpotent Lie algebra, 
denoted by §4, which is also a central extension of the Lie algebra Ej(3) by 
M. Indeed, taking the basis {ai, a 2 , a 3 , a 4 } of 04 such that the non-vanishing 
defining relations are 

[a x , a 2 ] = a 3 , [ax, a 3 ] = a 4 , (28) 

then the center 3 of the algebra is generated by {04}, and the factor Lie 
algebra 04/3 is isomorphic to ()(3). However, this extension is not equivalent 
to the extension appearing in the case of the planar rigid body with two 
oscillators, compare (|2^j) with (JTHJ). 

We give now briefly some results concerning the solution of (J25|) by the 
Wei-Norman method and its reduction to a system of Brockett type plus a 
system in R. The calculations are similar to that of the previous subsection. 

Let G4 be the connected and simply connected nilpotent Lie group whose 
Lie algebra is g 4 . The right-invariant system of type (0) on G4 corresponding 
to (j23|) reads as (JHJ) but where now g(t) is the solution curve in (5 4 starting 
from the identity, and {a%, a 2 , a 3 , a 4 } is the previous basis of g 4 . Writ- 
ing g{t) = exp(— vidi) exp(— v 2 a 2 ) exp(— v 3 a 3 ) exp(— ^404), the Wei-Norman 
method gives the system 

v 2 

vi = b 1 (t), v 2 = b 2 (t), v 3 = b 2 (t) Vl , v 4 = b 2 (t)^-, (29) 

with initial conditions fi(0) = v 2 (0) = v 3 (0) = v^O) = 0, which is eas- 
ily integrable by quadratures. The action $ of (5 4 on M 4 corresponding 
to the infinitesimal generators {X{\ reads, parametrizing g £ G 4 by g = 
exp(aai) exp(ki2) exp(cas) exp(da 4 ), 

$ : G A x M 4 — > M 4 
((a, b, c, d), fa, x 2 , x 3 , ar 4 )) 1 — > {xx, x 2 , x 3 , x 4 ) , 
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where xi = Xi — a, x 2 = x 2 — b, x 3 = x 3 — ax 2 + ab + c and x 4 = x 4 — ax 3 + 
a 2 x 2 /2 — a 2 b/2 — ac — d. The composition law reads 

(a, 6, c, d)(a', 6', c', <f ) = (a + a', 6 + b' , c + c' - ba', d + d' - ca! + 6a /2 /2) , 

the neutral element being represented by (0, 0, 0, 0). 
Thus, the general solution of (|23j) is 

$((-t>i, -f 2 , -v 3 , -V4), (xio, x 20 , £30, x 40 )) = (xw + «i, £20 + ^2, 

^30 + ^20 + ^2 - ^3, ^40 + ^30 + v\x W /2 + vfv 2 /2 - + V 4 ) , 

where t>i = t>i(t), t>2 = t> 2 (t), t>3 = v 3 (t) and t> 4 = f 4 (t) are the solution of 
(f2T?J) and the initial conditions are (x w , x 20 , x 30 , x 40 ) G M 4 . 

Due to the Lie algebra structure of g 4 , by quotienting by the center, we 
can reduce the solution of the equation in the group (5 4 (and hence of (J23)) 
to two other problems: one, a system in H(3) which is of Brockett type, and 
then we have to integrate a system in M. Taking canonical coordinates of 
first kind, and following analogous steps to that of the previous subsection, 
we have to solve first (j2*Uj) . and then 

d = bj f- (lyi(t)y 2 (t) - y 3 (t)\ - ^b 2 (t)yf(t) , (30) 
which is integrable by one quadrature. 
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